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In order to make plausible the idea that light exerts a pressure on matter, 

some introductory physics texts consider the force exerted by an electromag- 
netic wave on an electron. The argument as presented is both mathematically 
incorrect and has several serious conceptual difficulties without obvious resolu- 
tion at the classical, yet alone introductory, level. We discuss these difficulties 
and propose an alternative demonstration. 
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1 The Freshman Argument 

a: 

The interaction of light and matter plays a central role, not only in physics itself, 
but in any freshman electricity and magnetism course. To develop this topic most 
courses introduce the Lorentz force law, which gives the electromagnetic force act- 
ing on a charged particle, and later devote some discussion to Maxwell's equations. 
Students are then persuaded that Maxwell's equations admit wave solutions which 
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travel at the speed of light, thus establishing the connection between light and elec- 
tromagnetic waves. At this point we unequivocally state that electromagnetic waves 
carry momentum in the direction of propagation via the Poynting flux and that light 
therefore exerts a radiation pressure on matter. This is not controversial: Maxwell 
himself in his Treatise on Electricity and Magnetism[l] recognized that light should 
manifest a radiation pressure, but his demonstration is not immediately transparent 
to modern readers. 

At least two contemporary texts, the Berkeley Physics Course [2] and Tipler and 
Mosca's Physics for Scientists and Engineers[3], attempt to make the assertion that 
light carries momentum more plausible by explicitly calculating the Lorentz force 
exerted by an electromagnetic wave on an electron. In doing so the authors claim — 
with differing degrees of rigor — to show that a light wave indeed exerts an average 
force on the electron in the direction of propagation. Tipler and Mosca, for example, 
are then able to derive an expression for the radiation pressure produced by a light 
wave. 

A cursory look at the "freshman argument," however, which many instructors also 
present to their classes, shows that in several obvious respects it is simply incorrect 
and that in other respects it leads rapidly into deep waters. Nevertheless, one can 
more plausibly demonstrate that light exerts a radiation pressure and calculate it in 
a way that should be accessible to first-year students. It is the purpose of this note 
to discuss these matters. 

Consider, then, the situation shown in figure 1. We assume that a light wave 
propagates in the +z-direction, that its E'-field oscillates in the x-direction and that 
its Z?-field oscillates in the y-direction. The wave impinges on a stationary particle 
with charge q, exerting on it a force according to the Lorentz force law. In units with 
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Figure 1 : An electromagnetic wave traveling in the z-direction strikes a point particle with charge 
q. The -E-field is taken in the x-direction and the B-field is taken in the y-direction. 



c = 1, the Lorentz force is 



F = g(E + v x B) 



:i.i: 



which becomes 

F = q{E x - v z B y )\ + qv x Byk. (1.2) 

The "freshman argument" goes like this: Assume that E ~ sin(oot) and B ~ 
sin(uit). The particle is initially accelerated by the .E-field in the +x-direction and 
acquires a velocity v x > 0. The magnetic field then exerts a force on the charge 
equal to gv x B, which points in the +z-direction, the direction of propagation of 
the wave. The electromagnetic wave therefore carries a momentum in this direction. 
The Berkeley Physics Course in fact states, "...the motion of the charge is mainly 
due to E. Thus v is along E and reverses direction at the same rate that E reverses 
direction. But B reverses whenever E reverses. Thus v x B always has the same 



sign." 1 

A moment's reflection, however, shows that the last assertion is simply false. After 
one-half cycle, both E and B change sign. But because during this time the .E-field 
has accelerated the charge entirely in the +x-direction the electron at that point still 
has a positive x-velocity (In other words, the velocity and acceleration are 90° out 
of phase, as in a harmonic oscillator.) A similar argument holds for the z- velocity. 
Thus the cross product v x B reverses sign and now points in the negative z-direction, 
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opposite the direction of propagation. Furthermore, because there is an x-component 
to the force, one needs to argue that on average it is zero. 

The Berkeley authors indeed claim that the first two terms in Eq. (1.2) average to 
zero, the first because E varies sinusoidally, the second because B varies sinusoidally 
as well and because one "can assume that the increment of velocity along z during 
one cycle is negligible, i.e., we can take the slowly increasing velocity v z to be constant 
during one cycle." 2 With these assumptions the Berkeley authors conclude that the 
average force on the charge is (F) = q{v x B y )k. Although at first glance the result 
may seem plausible, it is also incorrect because the velocity and magnetic field are 
functionally orthogonal and consequently the time average of their product vanishes. 
That this is so, as well as the previous claim, can be seen by a proper integration of 
the equation of motion (1.2). 

2 Equation of Motion 

To determine the momentum of the charge, which we take to be an electron, assume 
the electric and magnetic fields of the light wave are given by E = E sin(uut + 0)i 
and B = B sin(ut + 0)j, where is an arbitrary phase angle. In our units E = B . 
Setting Fi orentz = mdv/dt in Eq. (1.2) then gives a pair of coupled ordinary linear 
first-order equations for the electron velocity: 

dv z 

— = u c sm(u)t + <p)v x 

— j- = u c sm(wt + <j>)[l - v z ], (2.1) 

where we have let qB /m = u c , the cyclotron frequency. 

These equations have the somewhat surprising analytic solutions 

v z (t) = c x cos — cos(ut + <fi) +c 2 sin 



U) c . , x 

— cosyut + <p) 

UJ 



+ 1, 
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v x (t) 



C\ sin 



UJ,. 



UJ 



COs{u)t + 0) 



c 2 COS 



UJ,, 



UJ 



COs{u)t + 0) 



(2.2) 



where C\ and c 2 are the integration constants. 

If we take 1^(0) = v x (0) = 0, which is reasonable and of sufficient generality for 
our purposes, we find 



C\ 



cos 



UJ r 



COS( 



UJ 



C2 



sin 



UJ r 



COSi 



CO 



and the full solutions are therefore 



v z (t) = — cos ( — cos 4> ) c °s 
v x (t) = — cos I — cos <p ) sin 



— cosf cut + 
a; 



UJ r , 



UJ 



cos(a;t + 0) 



— sin ( — cos 4> ) si n 

+ sin [ — COS0 ) cos 



UJ. 



— cosiujt + 

UJ 



+ 1, 



UJ r 



UJ 



cos(ujt + 



(2.3) 



The behavior of these solutions is not exceptionally transparent, but can easily be 
plotted. Figures 2-4 show several graphs for various values of uj c /uj and phase angle 
0. Notice that regardless of 0, v z is always positive, but that there is also a nonzero 
v x whose average can be positive, negative or zero depending on 0. The first case is 
shown in figure 2 and the third in figure 3. Also, in general v x » v z . 




Figure 2: The x-and z-velocities vs. cut for ui c /u> = .1 and = 0. Note that from the small io c jio 
approximation (see text), v x » v z always; in this case v z /v x = .1 
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^r- cot 



Figure 3: The same plot as in figure 2 except that 
average of v x = 0. 



= 7r/2. Note that in this case the time 




Figure 4: The same as plot 1 except that lu c /u> = .5 



Additional insight into the solutions can be obtained by examining the limit 
uj c /uj << 1. For ordinary light sources at optical frequencies oo ~ 10 16 rad s -1 , 
consideration of the Poynting flux (below) gives uo c /uu ~ 10~ n , and so the limit is 
well satisfied. For high-power lasers such as that at the National Ignition Facility, 
with energy ~ 2 MJ, one can have uo c > u. This limit should therefore be avoided. 
Expanding the solutions (2.3) to lowest order in uj c /uj for = yields 



i(£)W)-ir< 



— ) [1-cosM]. 



(2.4) 



Notice that both v z and v x are positive definite, as shown in figure 2. Therefore 
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their averages must be as well. This already contradicts the arguments of [2] that 
(F x ) = but that (F z ) ^ 0. Notice also that v z is of order (uj c /uj) 2 , while v x is of 
order u c /uo. This behavior coincides perfectly with the plots, but does suggest that 
since v x ~ (uj c /uj) 2 , a consistent, relativistic calculation will significantly change v z [4\. 
Moreover, the time averages of both v x and v z vanish to all orders, and so it is in fact 
impossible to exert a net force on the particle! 

One might object to the arguments of this section on the grounds that we have 
taken E and B to be simple harmonic ~ sin(u;£) rather than wavelike ~ sin(A;,2 — out). 
However, it is evident from Eqs. (2.4) that kz « uot always and that such corrections 
are therefore negligible, an assertion borne out by numerical calculations. 



3 Interpretation 

The question is now whether the behavior just discussed can be reconciled with the 
classical picture of the Poynting flux. The Poynting vector in our units is 

„ Ex B 



(3-1) 



4vr ' 

and the time average is (S) = Re(E x B*)/87r. S points in the direction of propaga- 
tion of the electromagnetic wave and in units with c = 1 can be regarded interchange- 
ably as power per unit area, energy per unit volume (or pressure) or momentum flux. 
If the freshman argument is correct, then the particle should be accelerated in the 
direction of the Ponyting vector, but our previous results show clearly that, to the 
contrary, the particle drifts off in some other direction at a constant average velocity 
and one searches for a way to explain away this fact. 

Unfortunately, there seem to be several deep inconsistencies in the entire ap- 
proach. A first is that the freshman derivation is evidently an invalid attempt to 
apply the standard classical derivation invoked to identify the Poynting flux with 
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electromagnetic momentum, a derivation which breaks down in the limit we have 
been considering. That is, advanced texts such as Jackson[5], typically begin by 
considering the Lorentz force on a volume of charges: 



dp 

~dt 



f (pE + J x B)d 3 x (3.2) 

J vol 



The first step is to eliminate the charge density p in favor of E via via Gauss's law, 
p = l/47rV • E. One also eliminates J in favor of V x B via Ampere's law to find 



"Pmecft , U J I f.^ ^ -d\ j3 



+ — / — (E x B)d s x 

dt dt J vol 4tt 

= -!- / [E(V • E) + E x (V x E) + B(V • B) - B x (V x B)}d 3 x. 

A.TT J vol 

(3.3) 

This is a purely formal result, which after the elimination of p and J relies only 
on vector identities. Since the second term on the left is the only one with a time 
derivative, one tentatively identifies it with the momentum of the field. 

However, the crucial difference between the "graduate" approach and the fresh- 
man method is that in the graduate approach we are considering a continuous charge 
distribution. In the limit of a single charge, the p in the Lorentz force law becomes 
the test charge distribution, whereas the p in Gauss's law becomes the source charge 
distribution and they cannot be equated. In the present situation there is not only 
a single test charge but no source charges whatsoever. Thus the standard derivation 
simply cannot be be applied. Indeed, the only volume one has at one's disposal is the 
volume of the electron itself, which leads quickly into quantum territory. 

A second difficulty is that the assumption of plane waves with constant amplitude 
is an assumption of constant energy and momentum. If the light wave has constant 
momentum, how can any be transferred to the electron? There are many instances 
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in physics where we ignore the backreaction of a recoiling particle on the system. For 
instance, according to conservation of momentum, a ball should not bounce off a wall, 
until one realizes that the ball's change in momentum is absorbed by the earth. 

Nevertheless, while to hold the amplitude constant in the current calculation might 
seem a reasonable approximation, to be totally consistent one should take into con- 
sideration the fact that the electron is accelerating and consequently emits radiation, 
and with that radiation momentum. The customary way to do this in the nonrela- 
tivistic limit is via Thomson scattering, but the differential Thomson scattering cross 
section for a wave polarized in the x-direction is 

-lk = \\ — I (cos 2 #cos 2 + sm 2 0), (3.4) 

d\l 2 \mj 

where 8 is the angle between the incident and scattered wave. 

The differential scattering cross section is defined as the ratio of the radiated 
power per unit solid angle to the incident power per unit area. We see that the 
Thomson cross section is absolutely symmetric with respect to reflection through the 
origin and consequently as much momentum is emitted in the forward as backwards 
direction. It is therefore far from obvious whether this situation can be corrected in 
the classical limit. Indeed, only when one goes to a quantum mechanical derivation 
(Compton scattering) does one see an asymmetry in the scattering cross section. In 
our situation, however, hu/m e ~ 10~ 5 , so it would appear that quantum corrections 
should be unnecessary. 

What one does practically to get the radiation pressure of light in, say, astro- 
physical calculations is to multiply the time-averaged Poynting flux (S) by the total 
Thomson cross section a^- One can see why this works as follows. A photon scattered 
off an electron will have a z-momentum p z = p cos 6, for initial momentum p , and it 
therefore removes (1 — cos 8)p from the original momentum component; the electron 
must gain the same amount. Multiplying the differential Thomson scattering cross 
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section (3.4) by (1 — cos#) and integrating over the sphere, gives exactly the total 
Thomson scattering cross section 

°T= B i( e -) 2 - (3.5) 

Multiplication by the momentum flux of photons will then give the total force on 
the electron. Because the Poynting flux is the momentum flux of photons the same 
numerical result is obtained by multiplying the Thomson cross section by the time- 
averaged Poynting flux. This entire argument, however, relies on the quantum nature 
of photons. The Thomson cross section is in actuality the nonrelativistic limit of a 
cross section that must ultimately be derived from QED, and so we see that the fresh- 
man plausibility argument leads quickly to a situation which may have no resolution 
in the realm of classical physics! 

The failing of Thomson scattering is due to the fact that no energy is removed from 
the original beam. A possible classical "out" to this situation is simply to assert that 
the energy radiated by the electron must be that lost by the incoming beam. There- 
fore since E = p for a classical wave, by conservation of momentum the electron must 
acquire a z- momentum exactly as in the Compton scattering case above [6]. While 
this argument may be valid in terms of conservations laws, it gives no mechanism for 
transferring the energy from the incident wave to the electron. Unfortunately, mod- 
elling the process as interference between the incident plane wave and the spherical 
wave outgoing from the electron fails to result in any transfer of ^-momentum from 
the wave to the charge. To recover the Compton result eventually requires including 
the radiation-reaction force on the electron, to which we now turn, but because this 
involves the classical radius of the electron it has already gone beyond the realm of 
classical electromagnetism. 

The most "straightforward" approach to deal with failure of the classical ap- 
proaches is via the Abraham-Lorentz model, which accounts for the energy radiated 
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by the electron, if in a somewhat ad hoc manner. From the Larmor formula the en- 
ergy radiated by an accelerated electron over a time T is ~ 2e 2 a 2 T/3. Equating this 
to the kinetic energy lost by the particle ~ ma 2 T 2 , one gets a characteristic time to 
lose all the energy to radiation: 

3m 
This timescale is 2/3 the time for light to cross the classical radius of the electron, 
r c = e 2 /m, and has a value r ~ 10~ 23 s. The total force acting on a particle will now 
be mv = F ext + F ra d, where F ra d is termed the radiation-reaction force. Conservation 
of energy considerations led Abraham and Lorentz to propose that F ra d = mr\ (see 
[5] for more details) and consequently the famous formula 

m (v - rv) = F ext . (3.6) 

With sufficient massaging, this equation can be applied to the present circum- 
stance to get the desired answer, that the force imparted to the electron by an elec- 
tromagnetic wave is F = (S)<tt- Eqs. (1.2) now become 

e 
Vx - rv x = —{E x - v z B y ) 
m 

v z - tv z = —v x B y . (3.7) 

In the nonrelativistic regime v z « 1 and we ignore the second term on the right 
in the top equation. For simplicity we also take both v x and v z to be of the form 
v = f e" lwt , which is of course manifestly untrue according to the results of §2. Then 
i) x = —iuv x and v x = —u 2 v x . The first of Eqs. (3.7) becomes 



J x 



or with lot << 1 



—iuivJl + icur) = — E x , (3.1 

m 



op 

v x = E x {l-iur). (3.9) 

muj 
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With the assumption that uj c /uj << 1 and lot « 1 we can ignore the v z term in 
the second of Eqs. (3.7). Then 

ic 
v z 2±——E x B y (l-iu,T). (3.10) 

For simplicity, take E x , B y real. Then we want the time average of the real part of 
this expression, or 

(F z ) = (mi,) = — 2 ^ = (S)a T , (3.11) 

as fervently desired. The earliest paper we have found that proposes this calculation 
is by Page [7] in 1920, although one suspects that Eddington carried it out earlier. 

Clearly there are a few things left to be desired in this derivation, but it does serve 
to show that the radiation-reaction force is necessary to get the claimed result. With 
slightly more work the conclusion can be put on a firmer footing via a perturbation 
calculation [6]: Note that Eqs. (2.4) are the zeroeth-order solutions of Eqs. (3.7), 
that is, when r = and v z << 1 is neglected. Assume v x = v x o + v x \ and v z = 
(v z q + v z i) « v x , where the subscript refers to the zeroeth-order solution and the 
subscript 1 refers to the perturbation. With the help of Eqs. (2.4)and (3.7) it is not 
too difficult to show that 

izi — u c v x i sin(ut) = uj^t sin 2 (cj£). (3-12) 

Taking the time average of this expression vindicates the previous result. We empha- 
size, however, that the Abraham-Lorentz model includes an explicit statement about 
the structure of the electron and hence cannot be regarded as entirely classical; the 
model is in fact a transition to quantum mechanics and quantum field theory. 

4 Alternative Approach 

Despite the many pitfalls revealed by the above methods, there is a superior and 
convincing demonstration that light exerts a pressure on matter, one that should be 
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accessible to freshmen who have had a basic exposure to Maxwell's equations. The 
great advantage of the method is that it avoids consideration of the force acting on 
a point charge and can therefore be carried out at the purely classical level. For 
this reason it should be adopted by introductory textbook authors. What follows 
is a simplified version of a calculation described by Planck in his Theory of Heat 
Radiation[8]. 




Figure 5: A light wave traveling in the z-direction strikes an almost perfectly conducting mirror 
of thickness z, width dy and height dx. An Amperian loop in the yz plane is also shown, with the 
direction of B given by the right-hand rule. 



Consider a light wave propagating, as before, in the +z-direction, and which 
bounces off a mirror, located at z = (see figure 5). We take the mirror to be a near 
perfect conductor of height dx, width dy and thickness z. The electric field of the 
light will now be a superposition of right- and left-traveling waves: 



E x = E cos(kz — out) — E cos(kz + out), 



(4-1) 



where k = 2tt/X is the wave number and where we have included a phase change on 
reflection. (This solution ensures that E = at the surface of the conductor. Recall 
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that the tangential component of an E-fi.eld must be continuous across a boundary, 
and since in the case of a good conductor the interior field essentially vanishes, the 
exterior field at the boundary must also.) 

From the differential form of Faraday's law, 3 V x E = —dH/dt we have 

V x E = — -] = -E k[sin(kz - out) - sm(kz + ut)]j = — — . (4.2) 

Integrating with respect to t and remembering that k = 10 in units where c = 1 
gives 

B = E [cos(kz — out) + cos(kz + ujt)]j = 1B cos(kz) cos(u;£)j. (4-3) 

Notice that at the boundary, B = 2B cos(ut) 7^ and that therefore by Ampere's 
law, § B • ds = 47r/, oscillating currents must be induced near the surface of the 
mirror. Since B is in the indirection, the right-hand-rule tells us that these currents 
will be in the ±x-direction, but that I x B will always point in the indirection. 
Therefore the Lorentz force due to the light, F = Idx x B for a mirror of height dx 
and total current / will in fact produce a force in the direction of propagation. 

We can calculate the magnitude of the force simply and plausibly. The magnitude 
of the Lorentz force is dF = IdxB, or dF = JdxdydzB, for current density J. Now, 
the differential form of Ampere's law tells us 

an 

V x B = -^p-\ = 4ttJ, (4.4) 

oz 

or J = —(l/A-K)dB y /dz. The Lorentz force therefore becomes 

dF 1 dB y 

-Bydz. (4.5) 



dxdy Air dz 



3 Most introductory texts use the integral form of Maxwell's equations. The derivation can easily 
be carried out by considering infinitesimal loops in the xz and yz planes as follows: The integral form 
of Faraday's law is § E • ds = —d<j>/dt for magnetic flux <p. In the case of our mirror, the right-hand 
rule gives 

/dW dT-i 

E • ds = E x (z + dz)dx — E x (z)dx = dzdx = — —dxdz, 

which leads immediately to Eq. (4.2). 

Similarly the integral form of Ampere's law § B • ds = 4nl leads to Eq. (4.4). 
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The quantity on the left is of course dp, for pressure p. Since the only spatial de- 
pendence of B is on z we can ignore the distinction between the partial and full 
differentials. Evidently, since dB y /dz is connected to J, we must interpret B as 
being the field exerting a force on a given slice within the conductor. Then, if we 
assume that the magnetic field drops off to zero at infinity, which is certainly true 
inside a good conductor where the falloff is exponential, the total pressure on the 
mirror should be 

P=--\ BdB = +-B(0) 2 = -Blcos 2 (ujt), (4.6) 

47T JO 07T 2ir 

where the last equality follows from Eq. (4.3) and the continuity of the tangential 
component of B across the boundary. The time average of this expression gives 

P — —, — 2(5') incident (4.7) 

47T 

as desired and where the factor of two is expected due to the recoil of the wave off 
the mirror. 

There are a few tacit assumptions in this derivation that should perhaps be made 
explicit. One might wonder, for example, why we use Ampere's law (4.4) to cal- 
culate the conduction current, rather than Faraday's law d<p/dt = — /E- ds = £, 
for magnetic flux <fi = Bdxdz and induced EMF £ . Normally, we would have 
students use this law to calculate the induced current / = £/R in, for example, 
a wire loop of resistance R. However, in a good conductor E « B and hence 
\d(f)/dt\ = § E ■ ds << /B • ds = 47rJ, the last equality representing Ampere's law. 
Furthermore, the 5-field in Eq. (4.4) includes both the incident field and that gener- 
ated by the induced currents. It seems unreasonable that the portion of the 73-field 
generated by the induced currents can result in a net force on the currents them- 
selves (no "Munchausen effect"). In fact this is the case and a detailed calculation 
demonstrates that integrated force on the induced currents by the induced 5-field 
vanishes. 
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Nevertheless, with these assumptions the simpler derivation we have presented 
appears sound, and it unequivocally shows that light waves do exert a pressure on 
matter in the direction of propagation. 

In conclusion we might say that, although one does not, and cannot, expect deriva- 
tions at the freshman level to be uniformly rigorous, this case is of particular interest 
because the interaction of light with matter is of fundamental importance. Moreover, 
the explanations presented in textbooks and in the classroom are so seriously flawed 
that even students sometimes notice the difficulties. Rather than try to paper over 
these problems with what must be regarded as nonsensical arguments, the occasion 
would be better exploited to point out that physics is composed of a collection of 
models that are brought to bear in explaining physical phenomena, but that these 
models have limited domains of applicability and, as often as not, are inconsistent. 
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